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Abstract 

A field theory describing the low-energy, long-wavelength sector of an incom- 
mensurate, spiral magnetic phase is derived from a spin-fermion model that 
is commonly used as a microscopic model for high-temperature superconduc- 
tors. After integrating out the fermions in a path-integral representation, 
a gradient expansion of the fermionic determinant is performed. This leads 
to an 0(3)(X'0(2)-symmetric quantum nonlinear a model, where the doping 
dependence is explicitly given by generalized fermionic susceptibilities which 
enter into the coupling constants of the a model and contain the fermionic 
band-structure that results from the spiral background. A stability condition 
of the field theory self-consistently determines the spiral wavevector as a func- 
tion of the doping concentration. Furthermore, terms of topological nature 
like the ^-vacuum term in (l+l)-dimensional nonlinear a models are obtained 
for the plane of the spiral. 

PACS numbers: ll.15.Tk, 71.27.+a, 75.10.Jm 
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I. INTRODUCTION 



The properties of doped magnetic systems have been intensively studied in the context 
of strongly correlated electronic systems like heavy-fermion compounds and cuprate super- 
conductors. The coexistence of local magnetic moments and itinerant fermions in these 
materials has raised many questions concerning the interplay between magnetic and charge 
degrees of freedom close to a magnetic instability. A powerful tool for a quantitative inves- 
tigation of the critical behaviour of given microscopic models is provided by renormalization 
group studies of the corresponding continuum theories. 

The field-theoretical approach has led to a successful description of the two-dimensional 
(2D) Heisenberg quantum antiferromagnet, where even a quantitative comparison with ex- 
perimental results for the undoped parent compounds of the cuprates was achieved Re- 
sults obtained from studies of the appropriate continuum theory, the 0(3) (2+l)-dimensional 
quantum nonlinear a model, through renormalization group analysis fl]] and chiral pertur- 
bation theory 0, agree very well with experimental 0] and numerical data. 

More recently, field theories for frustrated 2D Heisenberg quantum antiferromagnets like 
the one on a triangular lattice were derived. Since the starting point is the corre- 

sponding Neel state with a noncollinear order, the field theory is a (2+l)-dimensional S0(3) 
quantum nonlinear a model, as is generally expected for spiral states. Renormalization 
group analyses were performed for this model as well as for the classical 0(3) (8)0 (2) 
nonlinear a model in (2 + e) dimensions [|l^,|ri[]. Furthermore, the scaling properties of the 
quantum phase transition from the ordered helical to the disordered phase were described 



12|, based on scaling arguments and on a scenario of deconfined spinous. An understanding 



of the critical properties of frustrated spin systems between two and four dimensions was also 
achieved recently |13|. Thus, a fairly large amount of theoretical predictions was obtained 
in recent years in the case of pure quantum spin systems. 

However, the most interesting and relevant situation of strongly correlated systems, 
where magnetic as well as charge degrees of freedom interact, was until now investigated to a 
much lesser extent. Apart from numerous mean-field attempts to analyze those systems, to 
the best knowledge of the authors only few field-theories have been derived from microscopic 
models so far [p!^-[T7|, besides phenomenological approaches |]18[, where fluctuations effects 
are duly taken into account. 

The studies above were performed for holes in an antiferromagnetic background, such 
that no equivalent descriptions are available yet for the region of the phase diagram of 
strongly correlated systems, like some of the cuprates, in which incommensurate spin fluctu- 
ations appear. For doping concentrations in the superconducting regime, neutron scattering 



experiments have shown 1 19,20] that the Lanthanum compounds exhibit peaks in the mag- 
netic scattering intensity at wave vectors which are shifted by ±5(vr, 0) and ±5(0, tt) from 
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the antiferromagnetic point, where the shift 5 increases with the doping. The correlation 



length of the incommensurate fluctuations has been found to be considerably larger |20 



than the lattice spacing, which makes a continuum approach reasonable. Our aim here is 
to establish a mathematically rigorous connection between a suitable, sufficiently general 
microscopic model and a continuum theory. Considering fluctuations around a helical spin 
conflguration, we will develop a fleld theory that, in conjunction with the above-mentioned 
recent progress in renormalization group studies, leads to a quantitative description of the 
critical behaviour of an incommensurate spin system interacting with doped charge carriers. 

Previous theoretical studies of microscopic models of strongly correlated fermions suggest 
that the competition of hopping and exchange effects may lead to the formation of a spiral 
phase in the spin subsystem. Mean-fleld calculations favour spiral against antiferromagnetic 
configurations in Hubbard models for small but non-zero doping p^ , |2T[] . Concentrating on 
the one-hole case, Shraiman and Siggia proposed spiral spin ordering based on a phenomeno- 
logical continuum description of the t-J model |2^. They argued that the introduction 
of a low density of doped holes into a locally antiferromagnetic background may lead to 
an incommensurate, helical rotation of the staggered magnetization as a consequence of a 
coupling between the spin current of the holes and the magnetization current of the back- 
ground. Such a coupling was in fact derived in a continuum action obtained directly from 



the spin-fermion model IT^. Later on, the phase diagram of the Shraiman-Siggia model 
was investigated within a jj: expansion pSf. A transition from the commensurate to an 
incommensurate magnetic phase was obtained with increasing doping, both in the ordered 
and in the quantum-disordered regime for suitable values of the Shraiman-Siggia coupling. 
However, apart from the fact that the starting model is a phenomenological one, it is not 
clear to us to which extent those results really apply to the case of a spiral phase, since the 
conclusions on the critical behavior are derived on the basis of an order-parameter in the 
manifold S"^ and not SO (3), as it should be in the non-coUinear case. Finally, it should be 
mentioned that incommensurate spin configurations have further been shown to be induced 
by doping in numerical simulations of the one-band and the three-band Hubbard model 
away from half-filling [^J-|26|. 

The starting point of our work is the spin-fermion Hamiltonian [^, which describes mo- 
bile fermions interacting with a background of localized spins through an exchange term. It 
is a generalization of the Kondo lattice, since also a Heisenberg exchange interaction among 
the localized spins is included. We will concentrate on the two dimensional case that gives 
a realistic description of the cuprates. A central aspect that directs the choice towards the 
spin-fermion model is its being analytically tractable, an important advantage over other 
microscopic models for strongly correlated fermions. The study follows essentially the same 
steps as in the antiferromagnetic case [|15|,|T6|, although the actual calculation (Sec. IV) 
greatly differs from that case. In a path-integral representation of the model, the fermionic 
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degrees of freedom appear bilinearly and are integrated out exactly. An action containing 
only the spin degrees of freedom is obtained in terms of a fermion determinant and a pure 
spin part. We consider incommensurate, short-ranged spiral configurations for the spins 
and discuss their parametrization for long-wavelength, low-energy fluctuations around the 
ordered spiral phase. The gradient expansion of the fermion determinant is carried out in 
energy and momentum space in such a way that the occurring infinite series can be summed 
to all orders of the coupling constant by using the constraint on the order parameter. We 
show that in the limit of long wavelengths and low energies of the spins, the continuum 
theory is not only given by an SO (3) quantum nonlinear a model, as assumed in previous 
phenomenological approaches, but additional terms appear. On the one hand, a term linear 
in derivatives is obtained, such that parity is explicitly broken. However, since this term 
is not positive definite, it should vanish in order to guarantee a stable spiral configuration. 
This requirement leads to two equations that determine the wavevector for the spiral as a 
function of doping and the parameters of the model. On the other hand, a term of topological 
character like the one obtained in the continuum limit of the one-dimensional antiferromag- 
netic Heisenberg model p8| appears. However, in contrast to that case, the corresponding 
coupling constant is doping dependent, and hence, it can vary continuously. Therefore, our 
gradient expansion not only delivers the coupling constants of the SO (3) quantum nonlinear 
a model as functions of the microscopic parameters and generalized fermionic susceptibilities 
which contain the doping dependence, but shows important differences with phenomenologi- 
cal models obtained only on the basis of symmetry arguments. Finally, we would like to add 
that the present approach does not impose any restrictions on the energy and momentum 
scales of the fermions. 



II. MICROSCOPIC MODEL AND PATH-INTEGRAL DESCRIPTION 

We consider a spin-fermion Hamiltonian which describes spins localized on the vertices 
of a square lattice ("Cu-sites") interacting with fermions moving between sites situated on 
the bonds ("0-sites"). The problem remains nontrivial due to the coupling of the band 
fermions to the surrounding localized spins. The Hamiltonian is 

<jj'> i <jj'.i> 

a, a' 

<ii'> 

Here c^^ and Cj^a are creation and annihilation operators, respectively, for holes with spin 
projection a =|, | on sites situated on the bonds of the square lattice, denoted by the index j. 
The index i runs over the vertices of the square lattice. The kinetic term describes hopping 
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processes between the sites situated on the bonds, and it can contain a direct hopping 
between sites j and j' as well as an an effective hopping between sites j and j' mediated by 
a central site i. The term proportional to Jk is a non-local, Kondo-like interaction between 
the localized spins Si and the holes on the neighbouring sites. It includes spin-exchange 
processes with hopping in addition to pure exchange processes. The vector (Taa' consists of 
the three Pauli matrices. The summation denoted by < jj', i > runs over all pairs of sites j, 
j' which are nearest neighbours to a given site i. Finally, ifgf contains an antiferromagnetic 
Heisenberg superexchange interaction between nearest-neighbour spins Si. 

The spin-fermion model can be obtained as the strong-coupling limit of the three- 



band Hubbard model [^, which, according to numerical simulations p5| , |30| , consistently 
describes the Cu02 planes of the cuprates. Limiting cases of the spin-fermion model cor- 
respond to other models that are frequently discussed in the context of strongly correlated 
fermion systems. In the limit of zero doping, the spin-fermion model reduces to an antifer- 
romagnetic S=| Heisenberg model, which yields a quantitative description of the undoped 
cuprates [0,0,3 • For large Jk, the doped holes on the O sites strongly bind to the central 
Cu ion to form a local singlet, so that for zero direct 0-0 hopping the low-energy dynamics 
of Hgf can be mapped onto the t-J model ||3l|. For finite 0-0 hopping, the spin-fermion 
Hamiltonian can be mapped onto a generalized t-J model containing second and third 



nearest-neighbour hopping and spin- flip hopping . For vanishing Heisenberg interaction, 
Hsi is equivalent to a Kondo-lattice Hamiltonian with a nonlocal exchange between band 
fermions and localized spins. 

In the following we discuss only briefly the path-integral representation of the model in 
order to set up the notation, since the same steps were already performed in the antiferro- 
magnetic case \T^\. We represent the partition function of the spin-fermion Hamiltonian as 
an imaginary-time path-integral. Employing a spin-| coherent state representation |Q for 
the spin degrees of freedom and Grassmann variables for the hole degrees of freedom, we 
obtain for the partition function of Hgf 



^sf 



y ©5 j Vc*Vc exp(Ss + 5f) 



(2) 



S(0)=S(/3) c(0)=-c(/3) 

where 13 = l/ksT. In Eq. (|^), the pure spin part of the action is given by 



dr 



■J2A{S.{T)/S)-d^S.,{T)-JnJ2 



<ii'> 



(3) 



where now S = Sfl and fl = (sin^,;(r) cos</),;(r), sin0j(r) sin0j(r), cos^^j(r)). As is well 
known, the first term in the action Ss is the Berry phase for the adiabatic transport of a 
quantum spin along a closed circuit and responsible for the correct quantization of the spins. 
The monopole potential A satisfies the constraint 
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abc^ = fi« (4) 

on the unit sphere. The measure for the integration over the spin variables is given by 
In momentum and frequency space, and after diagonahzing the kinetic part of Hgf, which 



yields two bands with a dispersion relation e(fc, A), A = 1, 2, for which we refer to Ref. ||T5 
we obtain 

'S'f = /3 ^ ^ ^ cl,kx (i^n - e(fc, X)^6aa> 6kk' Sxy 

a,a' k,k' A,A' 

-4Jxcraa' ■ Sk-k'S*{k, A) s{k', A') Ca'k'X' , (5) 

where k = (e„,,fc). Its zeroth component is a discrete Matsubara frequency defined by 
e„ = {2n — 1)71 /f3 for fermionic fields, and by Un = 2mT/f3 for bosonic fields. The form 
factors s{k, A) contain information on the band structure of the free system and are given 
by s{k, A) = ei(fc. A) sin {kxa/2) + e2(fc. A) sin {kya/2), where ei^2{k, A) are the components 
of the eigenvectors ||T5[ that diagonalize the kinetic part of Hsf. 

Since the action in Eq. (|^) is bilinear in the Grassmann fields, the integration I)c*I)c 
may be carried out within the path-integral to obtain 

Zsf= J VSexp{S, + Ti\nG-^) , (6) 

S(0)=S(/3) 

where 

G-^ = Go 1 - S (7) 

is the inverse propagator of the fermions in the presence of the dynamical spin field Sk-k'- 
The free fermionic propagator is given by 

iGo^)akX,a'k'\' = i^^n " e(fc. A)) 5aa'5kk'SxX' , (8) 

and S is the self-energy of the holes interacting with the Cu spins, 

i^)akX,a'k'X' = (Taa'Sk-k' s\k, \)s{k', A') . (9) 

The trace in Eq. (^) is to be taken over the indices k, a and A. It is convenient to shift the 
form factors s{k, A) from the self-energy S to the unperturbed fermionic propagator Go- K 
can be checked that by writing the logarithm as a power series, Tr ln(Go ^ — S) = Trln(Go ^) — 
Yl'^=i n'^^(^o^)"' rearranging the terms in the matrix product, a redefinition of the free 
fermionic propagators 
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A — 1 



and the self-energy 



i'^)ak,a'k' - CT aa' ' Sk-k' , (11) 

yields the relation 

2 

Tr ln(Go^ - S) = Tr Infj^ |.(fc, A)rGo ^(3 - A)) + Tr ln(Go ^ - S) . (12) 



Since the first term on the right-hand side of Eq. ([T^) contains only the free propagator 
and no interaction contribution, it may be regarded as a normalization constant and will 
be ignored from now on. Note that the trace in the second term on the right-hand side of 
Eq. (0) no longer includes the band index A. Since in the following we will use the free 
fermionic propagator and the fermionic self-energy only in the forms (0) and ([lT|) , we will 
omit the bars to simplify the notation, and refer to Eqs. (p!OD and ([TT|) simply as Go and S. 
We further introduce the abbreviation 

^°(') = 5:(.,„_e(fc,A))- ^''^ 



A=l 



Later on we will need the momentum-space symmetry properties of the function go{k) = 
gQ{en,k). Using the explicit expressions [0 for the energy eigenvalues and eigenvectors of 



the free fermion system, the following symmetry relations can be readily deduced: 

go{en, k) = go{en, -k) , (14) 
go{en,K,ky) = go{en,-K,ky) = go{en, k^, -ky) . (15) 

We have derived an action where the only degrees of freedom appearing explicitly are the 
spins, which are suitable for a continuum approximation. This continuum approximation 
does not affect the fermionic degrees of freedom which are taken into account with their full 
dispersion relation. 



III. ORDER PARAMETER FOR THE SPIRAL 

For the spin configuration Si along the lattice, we have to introduce an expression char- 
acterizing the spiral order expected in the classical ground state. We can then take into 
account long-wavelength, low-energy fluctuations around this ordered helical state. In the 
classical ground state, the spins lie in a plane in spin space. In order to derive a continuum 
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theory, we need to identify vector fields which are smooth in the long-wavelength, low-energy 
regime. In the case of a Heisenberg antiferromagnet on a triangular lattice, the constituent 
fields of the order parameter can be obtained from linear combinations of the physical spins 
within one magnetic cell PJ^, in analogy to the construction of the order parameter fields 



describing coUinear antiferromagnet configurations |3J]. In the case of an incommensurate 
helix however, there are infinitely many sublattices and a magnetic cell does not exist. 

In order to construct the order parameter for the spiral, we consider long-wavelength 
fluctuations around a classical spiral configuration with wave-vector 

ks = Q + 6Q, (16) 

where 

«=(-■-) (IT) 



.a a . 

is the wave vector of the antiferromagnetic state, while 

6Q = i6Q.,5Qy) (18) 

is incommensurate with the lattice, i.e. m6Q ^ nQ for all m, n integer. We start with 
the vector product of neighbouring spins. We define first a vector n^i = Si x Sj+i /Ni, 
where Ni = S"^ sin0j_j+i and is the angle enclosed between the spins Si and Si+i. 

This choice is well defined as long as the line joining the two spins is not perpendicular 
to kg. The vector reduces in the classical ground state to a vector parallel to the 
axis of the helix, which is independent of the lattice site. We now define a matrix i?"^ = 
cos{ks-ri) 6°''' + [1 — cos(fcs-Tj)] n'^i'n^i — sin(fcs-rj) e^'^^rig- which performs rotations about 
risi. Here the upper indices a,b,c = 1, 2, 3 denote components in spin space. We adopt 
the convention that summation over repeated indices is implied, although for clarity we will 
occasionally write the summation sign explicitly. Applying Ri to the spin Si rotates it back 
into a space-fixed axis, which provides the definition for a second vector being constant in 
the classical ground state: riu = Ri Si / S. Defining a third vector n2i = n^i x Uu, we have 
a complete set of orthonormal vectors which are constant in the classical ground state and 
fulfill 

nlnl = 5^,. (19) 

Using the above basis, we parametrize the incommensurate spiral configuration as 

Spq{T) = S [ni{p, q, t) cos {k.-Tpq) - n2(p, g, r) sin (fc^-r^g)] , (20) 

where now p, q (instead of i) denote the sites on the two-dimensional square lattice and 
fpq = {p ' ' o,), with a the lattice constant. It is evident from Eq. (|20|) that in the 
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limit 5Q —>■ 0, the antiferromagnetic ground state is reproduced, so that 5Q describes the 
deviation from perfect Neel order. 

In order to evaluate the continuum limit of the action, we allowed the vectors rii and n2 
to be smooth functions of the lattice site and of the Euclidean time r. The constraint (|T9D 
must be satisfied at every given point in space-time. In a previous work [^, we have treated 
the special case of planar fluctuations, where the variation of ni{p, q, r) and n2(p, q, r) was 
confined to the plane which is defined by rii and n2 in the classical ground state. The order 
parameter for the description of planar fluctuations around a spiral configuration is a two- 
component unit vector in the manifold 5*1. The corresponding field theory was shown to be 
an 0(2) nonlinear a model. Here we extend this treatment to the general case of fully three- 
dimensional fluctuations. This generalization turns out to be highly non-trivial (see Sec. IV 
and V). In the present case, fluctuations of ni(p, g,r) and n2(p, g,r) imply fluctuations of 
n3(p, g,r), so that the excitations are described by an S0(3) order parameter: 

Q = (ni 71,2 ns) or Qab = < • (21) 

In contrast to an antiferromagnetic configuration, where the ground state is invariant under 
rotations about the spin axis, in a noncoUinear ground state global rotations about any axis 
lead to a new ground state. Thus in the spiral phase the global 0(3) rotational symmetry 
in spin space is completely broken and we expect three Goldstone modes. 

In helical phases, however, there is still a local rotational symmetry, which was first 
identified in Ref. S'pq(r) is invariant under rotations by an arbitrary local angle ippqir) 



about the /oca/ axis Spqir). This can be readily seen in a Schwinger boson representation of 
the spins where tu" is a doublet of complex scalar fields satisfying uj'^lu" = 1 

Il36[| . In this representation, a local rotation about the local spin axis corresponds to the 
U(l) gauge transformation However, under these gauge transformations the 

order parameter fields rii = 1, 2, 3) are not invariant. In the general case, in which 5Q 
is finite, an infinitesimal change in the angle ippqir) between two sites will lead to a finite 
change in the fields n^. If the fields rii are slowly varying in one choice of gauge, they will 
have rapid variations in other gauges. Thus, by focusing on continuous configurations, one 
explicitly breaks this local symmetry by fixing the gauge. Our prescription to determine the 
fields rii, and therefore to fix the gauge is physically natural, since it is given by the actual 
configurations of the spin fields S. 



IV. GRADIENT EXPANSION FOR THE DOPED SYSTEM 

In order to parametrize the low-lying modes around the ground state of the classical 
spiral, we decompose the spin field into a helical and a uniform component so that 
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2r2\-l/2 



S = S {n + aL) {1 + 2 an ■ L + a^L'^) 
where we have introduced the abbreviation 

n = ni cosiks ■ Vpg) - n2 sin{ks ■ Vpg) 



(22) 



(23) 



and the {p, q, r)-dependence of the fields is imphcit. X is a slowly varying ferromagnetic 
field with \aL\ <C 1, which corresponds to the net magnetization density. L is of the same 
order as a first-order derivative of n. While n describes fiuctuations around the spiral wave 
vector kg, L describes fiuctuations around k = 0. Since the resulting action will be bilinear 
in the ferromagnetic fiuctuation, L can be integrated out at the end of the calculation. The 
inverse square root factor in Eq. (^) puts S/S on the unit sphere. 

Expanding the ansatz for the spin field ( p2D up to second order in a gives 



S = S <n + a 



{n ■ L)n 



(n-L)L 



-L^n--(n-Lfn\} 
2 2^ ' \) 



(24) 



2 2 

This expression will now be employed for the space-time-dependent spin field in the path- 
integral, Eq. d^). We will perform the gradient expansion of the pure spin part of the action, 
Eq. (|), in real space. It is technically advantageous to carry out the gradient expansion 
of the fermion determinant (see Eqs. (|^)-(|3)) in Fourier space, where it is approximately 
diagonal for long-wavelength spin fields. We therefore transform the expression Eq. (|2^) for 
Spgi^r) into /c-space (where k = (e„,fc)) and insert it into the fermionic self-energy S (see 
Eq. (|l3))- This leads to an expansion of the fermionic self-energy in powers of a: 



Introducing the abbreviations: 

n± = n-j_ ± i n2 , 

we obtain for the zeroth, first and second orders of the fermionic self-energy: 



' ak.a'k' 



9 a 
2 '^aa' 



(E«) 



ak.a'k' 



rf_{k-k' -ks) + n\{k-k' + ks) 
9<a'[L''{k-k')-\Y,L\q,)x 

91:92 

[ 5^ <(g2 + rA;,)<(A;-fc'+rA;,-Ej^,g,) 

r=-,+ 

+ 2 5^ n\{q,)n'^,{k-k'-Y.U1^)]\^ 

d=l,2 ^ 

i^^"^)ak,a'k' =^?<«'{-^E^'^(^i)^'fe) E ^ik-k'+rkg-Eli^d 

^ „. „„ r=-,+ 



(25) 



(26) 



(27) 



(28) 



91.92 
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2 

d=l r=-,+ 
2 

d=l r=-,+ 

+ 5Z ri^(g3 + rM<(94 + rfc,)<(A;-fc'+rfc,-Eti9i)] k (29) 

r=-,+ ' 

where we have defined the couphng constant 

9 = 'iJKS. (30) 

Note that kg = (0, fc^). Inserting this expansion of the fermionic self-energy into the fermion 
determinant from Eq. (j^), it may be written as: 

Trln(G'o^ - S) = Tr ln(G'o ^) + Tr ln(l - aGoS^^^ - a^GoS^^)) , (31) 

where 

Go 1 = Go ' - . (32) 

The first term in Eq. ( |3TD represents the hehcal contribution since the field L does not enter 
in it. The second term contains the contributions of the ferromagnetic field. The helical 
and the ferromagnetic contributions to the fermion determinant, Eq. (|3T1) , will be evaluated 
using different methods, which both rely crucially on the spin-field momenta being small. 

A. Spiral contribution to the fermion determinant 

As a first step, the logarithm in the first term on the right-hand side of Eq. (^) is 
expressed as a power series: 

oo _ 

Trln(Go ' - S(°)) = Trln(Go ') - ^ -Tr(GoS(°))" . (33) 

The first term represents the free part and will enter the path integral as a multiplicative 
normalization constant. Our task is to take the trace over spin and momentum indices, and 
obtain a closed expression for general n, so that the infinite series can be re-summed. It is 
useful to write down explicitly the trace over the matrix product in the above equation: 
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Oil. ..an 



9o{kn)aZa^ [n^-ikn - h - k^) + - fcl + fc. 

Multiplying out the brackets leads to a sum 

J] ... ^ (7o(A;i)<,,<(^l-^2 + rifc.)...(7o(fcn)Cm<:(^n-^l+^nM- 



(34) 



(35) 



Since we are taking the trace, the sum of all momentum transfers to the spin fields must 
be zero, so that A;„+i = ki. Since kg = (0, ks) (where kg = Q + SQ, see Eqs. (H)-(|ll)) is 
incommensurate, this can be fulfilled only if n is even and the number of r j = (— ) equals the 
number of rj = (+). Consequently, the summmation ^ over all configurations of the 
{vi} contains (^"2) terms, namely all permutations of {n/2) n_-fields with {n/2) n+-fields. 
Since n is even, the trace in spin space over the string of Pauli matrices can be carried out 



using the trace reduction formula p7| : 

Tr {a^'a"^ ... a"") = 2 ^'(-l)^5"i"'2 . . . 



(36) 



(37) 



where P is the permutation 

and the sum Yl'p includes permutations between different index pairs only. In order to 
perform the gradient expansion of the expression Tr(GoS*^''))", the momenta in the arguments 
of the fields are redefined |^| in such a way that the free fermionic propagators go appearing 
in Eq. ( [3^ ) can be expanded in powers of the momentum transfer to the spin field. We obtain 

Tl' (G„S<»')" = 2 



2/ E EE' 

(Ji...g„_i ri...r„ P 



X 



^ goik)go{k -qi + nks) . . . go{k - + Y.i=ii^i ks) 

k 

<teK2fe-gi)...<:(-gn-i) 

We abbreviate the product of propagators as 

n(gi . . . Qn-i, ri . . . r„) = ^ go{k)go{k -qi + riks) . . . go{k - + YH=i^'i ks) 



(St 



(39) 
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Note that the form of the function 11 depends on the set of parameters {vi}. The g's are 
small, so we can expand the function 11 around k. It should be stressed that only the 
momenta exchanged with the spin fields are small but not their values for a given fermionic 
state. By Fourier transformation of the spin fields back into direct space and subsequent 
integration by parts, one obtains from an expansion up to 0{q^) first and second order 
spatial and temporal derivatives of the spin fields. Defining 

rf = /" d\dr ^ [nil (^) • • • dXr\ (^) • • • <: (^)] (^) • • • <: (^) > (40) 
1=1 

dXri i^)--- d.Kz (^) • • • < (^)] <:i (^) • • • <: (^) ' 

(41) 

= J d'xdr { ^ [Ki (x)... d,Kl (x)... d,K2 (^) • • • K: i^)] K:ti (^) • • • K: (x) 

l,m=l 
1=1 m=i+l 

we have 

Tr (GoS^"))" = 2 (-)" ^ 5^ J^'i-lfS"'"^' ■ ■ ■ ^'^'-i"- 
V 2 / ci 

ri...r„ P 

Here 9^ is an abbreviation for d/dx^. Analogously to our notation for the 3-momentum, 
in which k = {en,k), we define x = {t,x), so that the indices /i and u in Eqs. P0|)-(^3D 
run over three values. In order to avoid confusion with the summation index i, we use i to 
represent V— 1 in this section. The boundary terms resulting from the integration by parts, 
which we performed to arrive at Eq. (|4^), vanish because the field n{x) has been required 
to be constant in infinity, which is a natural assumption in a low-temperature approach. As 
is evident from Eq. (|43|) , the Kronecker deltas 5"'"-' generate a pairwise contraction of the 
fields and their derivatives into inner products. Carrying out the sum we sum over all 
possible pairwise contractions. In the derivation of Eq. (^31), the local part containing the 
zeroth order term in the expansion of 11 has been discarded since it contains no derivatives 



I A'xdr ^ [n'll 

Lm=l 



X) 
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and is therefore constant due to the constraint, Eq. (0). Several types of inner products 
are encountered in Eq. (|43D : 



nlnl = n''_n''_ = , (44) 

(9^<K = (9^n^K_ = 0, (45) 

< = 2 , (46) 

{d^nt_)n\ = -{d^nX)nt_ = 2 i {d^nl)nl , (47) 

{d,n-_) {d,n\) = {dXi) {d,nl) + {d,n'^) {d^nl) . (48) 

Pair contractions of the type {d^n"i) {di^nt) need not be considered since they imply a con- 
traction of the type n^n" (because the number of fields is even), which gives zero according 
to Eq. (|4^). From Eqs. (|4^)-(p8D, it is clear that one gets nonzero contributions only from 
those contractions 5""'""^ for which = — rj. Since there are (|)! possibilities for such con- 
tractions producing inner products of (|) n+-fields with (|) n_-fields, the sum in 
has (|)! terms. 

We will now outline the combinatorial analysis which we have performed in order to 
permit an explicit evaluation of the sums appearing in Eqs. (^ ) - (|3) to arbitrary order 



in n. This analysis departs considerably from previous work in the antiferromagnetic case 
T5| , p!6| . We denote the contribution of the first-order derivatives to Tr(GoS*'°'')" by 



T\...V- ^ 9—1 -^^ 



where was defined in Eq. ( |4D| ) . We will show that the only nonvanishing contributions to 
the sum in X!^ come from strictly alternating configurations of the set {rj}, i.e. those 

configurations which fulfill 

Let us consider a nonalternating configuration of the {rj}, where at least two ?24.-terms 
(and consequently two n_-terms) are nearest neighbours, i.e. we assume rj = rj+i for some 
i. One of the possible contractions connecting only ?T,+-fields with n_-fields is: 

[. . . n!! . . . n"^n"^+' . . . n!.^ . .] . (50) 

Another contraction connecting only n^-fields with n_-fields by is obtained by one permu- 
tation: 

[. . . n!! . . . n"^n"^+' . . . . .] S^^^j+iS^^^>^ . (51) 

Two sets of contractions differing by one permutation carry the sign (—1)-'^, (— l)^"*"^, re- 
spectively. Since the two terms ( ^OD and ([5T|) are identical except for their relative sign, 
they add up to zero in the sum Yl'p- Since (|)! is an even number, one can find for every 
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nonvanishing set of contractions another one differing only by a relative sign given a nonal- 
ternating configuration of the {^j}. Thus all terms in the sum add to zero pairwise for 
a nonalternating configuration of the {vi}. From the sum we are left only with the 

two alternating configurations. 

For a strictly alternating configuration of the {vi}, where rj+i = — rj V i, the total number 
of nonvanishing terms from the sum Y2'p is obtained by the following consideration: If a 
pair-contraction 5°'°^ with = —rj leaves a nonalternating sequence of the {vi}, all terms 
will again add to zero pairwise. The only contractions that always leave an alternating 
sequence of the {rj} are contractions between nearest neighbours or between the first and 
last member of the remaining sequence. The number of these contractions is 2^^~^\ and 
they carry the same sign. Thus in Eq. (^91) , the sum Y2'p{—1)^6"'^'^'2 . . . S^^"-'^"'^" can be 
replaced by 2^5~i) times the contribution of the zeroth permutation, §"-^"-^§^304 _ _ _ ^a„_ia„_ 

After the sums r„ Yl'p have been discussed, we still need to carry out the sums 
J2i=i Si=i (s^^ ^'i^- order to obtain a closed expression for X^. First we 

consider the alternating configuration starting with ri = (— ). For this configuration of the 
{rj} and the set of contractions §"'^"'^6°''^°"^ . . . ^'^^-i'^"^ every term with / odd gives {dnt)n°;^, 
while every term with / even gives {dn'^)n°i. There is always one contraction containing the 
derivative, for which we use Eq. p7|), and (| — 1) contractions without derivative, for which 
we use Eq. (|46|). This leads to 



(f-i) 

1=1 1=1 

1 for i odd 



2ti(90n^ X { ^ " (52) 
for i even. 

For i odd, the first derivatives of the product of propagators 11 (defined in Eq. (|39D) for an 
alternating configuration starting with ri = (— ) acquire the form 



l-odd 



dgo,~{k - q) 



g,;=0 ^ Oq 



g=0 



^EMfc)]^^k,-(fc-g)]%^^. (53) 



n ^ dqi^ 

k 



We define 



9o,±{k) = go{k ± K) . (54) 



Since dU ^ Idq't L =o has the same form for all i odd, the summation over i can be carried 

' ■'■'■odd I "i* ' 

out: YllZiiodd^ ~ f- Inserting these results into Eq. (^9]), and recalling that the sum Y^'p 
over the permutations is replaced by 2'^t~i) times the contribution of the zeroth permutation, 
we obtain 
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I n dli ^ 



n „9n- 



2a2^ dq^i 



odd 



qi=0 



q^=o 

d^xdr {d^nl)n^2 



(55) 



For our fermion determinant, Eq. (^), we need —Yl'^=i n-^n{-+)- insert Eq. (^3]) into 
Eq. (^) and re-sum the power series in n: 

n—9 n = l h 



n — 2 

n even 



n=l ~ k 

1 d 



9=0 



2a2 ^q^' 



^ln[l - g'go{k)go,.{k - q)] j ^^xdr {d^n\)nl . 



(56) 

Adding also the contribution X^^^^ from the alternating sequence starting with ri = (+) 
and introducing 

$f (g) = ln(l - g'go{k)go,±{k - q)) , (57) 



we obtain 



^ n 



2a? dq^ 



(<l>+(g) - $r(g)) j d^xdr (9^<)n^ . 



(5^ 



The sum in our expressions for ^"^{q) denotes a sum ^. By shifting the summation 
index k —k in ^f{q) and using the symmetry of the free propagator, Eq. (|14|), it can be 



shown that 



dqO 



q=0 



dq 



1,2 



dqO 

d^^iq) 



-,1,2 



q=0 
q=0 



(59) 



q=o dq 

such that the time component of the derivative drops out and only the space component 
remain: 



n=2 



q=0 



d^a;dr ni ■ 5^n2 . 



(60) 



For SQ = 0, i.e. a pure antiferromagnetic configuration, the first-order contributions vanish. 
We will discuss the relevance of this term for the resulting field theory at a later stage, after 
we have also obtained a contribution linear in the derivatives when we take the continuum 
limit of the Heisenberg part of the action. 
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We now turn to the second-order terms in the expression (^31). The contribution of the 
quadratic derivatives is denoted by 



i=i " 



(see Eq. (^) for F^^). We proceed similarly as in the case of the first-order derivatives, and 
investigate first whether for some special configurations of the {rj} the terms in the sum 
over all contractions Yl'p ^'^^ zero. It will turn out that in the case of the quadratic 
derivatives not only the strictly alternating configurations of the {r-j} give a nonvanishing 
contribution, but also configurations where the alternation is broken once. Again, in a 
nonalternating sequence one can find two possible sets of pair-contractions, differing by 
one permutation, that connect only (H — )-pairs and carry the signs (—1)^ and (— l)-^"*"^, 
respectively. The corresponding terms will only add to zero if they are identical except for 
their signs. However, this is not always the case if we deal with the second-order-derivatives. 
It may occur that in one of the two sets the two derivatives are contracted, 

[. . . a^n!: . . . d^n"^ . . . n+'= . . . n!! . . .] 6"'"^ 6'"''" . (62) 

while in the set differing by one permutation the derivatives are contracted with other fields, 
i.e. 

[. . . a^n!! . . . d^n"^ . . . r^+^ . . n!! . . .] S^'^'^'S''^'" . (63) 

While a contraction of type ( p^ ) gives a term of the form [(9^ni) ■ (duTii) + (9^n2) ■ {dun2)] 
(see Eq. (^8])), the contraction (|63D yields 4(ni ■d^n2){ni ■dyn2) (see Eq. (|47|)). These terms 
are not identical but related through 

2{ni-d^n2){ni-dyn2) = (5^ni) ■ (d^rii) + (S^na) ■ {d^n2) - (S^ng) ■ (S^ng) . (64) 

Whether two terms carrying the signs (—1)^ and (—1)^+^ are identical depends on the 
position of the derivatives, which is denoted by / and m (see Eq. (^)) and varied in the 
sum ^Y^i^^Y Thus we will now evaluate the sum for all possible positions / and m of 
the derivatives and then compute the sum ^^^^i- Again, we first discuss the alternating 
configuration starting with r\ = (— ). The contribution of terms with / = m is 

Ttiln.) = E' (-l)'^^"'^- • • • 5'^--"'" / d'xdr [Kl . . . (9,9.<;) . . . <:] 
p 

2 

= -2("-2) Y.^d,n,) ■ {d^n,) . (65) 

b=l 

For terms with / ^ m, where /, m are both even or both odd, we obtain (2;eN\{0}) 
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Terms with / = m ± 1 yield 



(66) 



(/=m±l) 



p 

2 

= 2("-3) ^(a^n,) ■ (9,n,) + 2("-=^) 2 (ni ■ d,n2){n, ■ d^n^) . 

b=l 

Terms with / 7^ m ± 1, where / even, m odd or vice versa, contribute 

p 

■ ■ ■ \^''"'ri±2z+i) ■ ■ ■ "'r„ 

2("-2) 2 {m ■ d^n2)in^ ■ d,n2) . 



(67) 



^ {l=m±2z+l) 



We now determine the number N of each of these four types of T^'^-terms in the sum 

Y!i,m=i^ denote the result of the summations [X)z,m=i Y^p] by i-e. % = N(^i=m) T^il^n) + 

iV(z=m±2.) 7^(z=^±2.) + N^i=m±i) ^^(Zmii) + N(i=m±2z+i) Tj^L^^iz+iy Using Eq. (0), we obtain 
that % does not depend on explicitly on namely. 



7;(z = odd) = -2(^^-2) ^(a^rife) • {d,n,) 



(69) 



For the alternating configuration starting with ri = (— ) which we have discussed so far, the 
second derivative of the product of propagators from Eq. (BUf) becomes 



dq'^ dqf 



1^ dq^ 



qi=0 



gi=o 



J2 9o{ky^~' 9oAk) 



„ 92^0 (^ - q) 



<7=0 



dq'^dq" 



q=0 



(70) 



so that the sum Y17=i '^^^ be carried out. Adding the contribution of the alternating 
configuration starting with ri = (+), we obtain for the contribution of the two alternating 
configurations to the sum ^ in Eq. (|6TD : 



II 

n{alt) 



2{n-2) 



9o{ky 9o,-{ky 



^id'^go,-{k- q) 
dq^^dq'^ 



q=0 
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--1 

2 9o, 



n d^go{k - q) 



dq'^dq" 

2 

d'^xdrl - ^((9;. rife) • (d^rii,) + 



q=0 



X 



n 
2 



1) {d^ns) ■ (a^ng) 



(71) 



For a nonalternating configuration, the terms Tj^^^^-^ aii^i '•-{i=rn±2z) 



vanish. The terms 



and T, 



and T|/^'' 

./, anKj. J.,, I o , IN can be nonvanishing when the contraction of the two derivatives is 
such that the remaining sequence is alternating, i.e. when the contraction of the derivatives 
'repairs' the defects in the alternation. However, such terms can still add to zero in the 
sum A term with fixed /, m, where the contraction of the two derivatives restores 

the perfect alternation of the remaining sequence, adds to zero e.g. with the term /, m + 1, 
where the 'repairing' contraction differs by one permutation and thus carries the opposite 
sign. Therefore, from both summations ^^^^^ and X]p) ^^^^ S^^ ^ nonzero contribution 
only if just one single term, where the contraction of the derivatives restores the alternation 
of the remaining sequence, occurs in This happens precisely if the alternation of 

the configuration is broken only once. These configurations will be denoted by 'one-kink 
configurations' in the following. 

We proceed with the discussion of the summations Yl^Zi-, Yl\m=i-: Y^p- denote 
the site after which the kink occurs by t. From the sum 'Y^Zl only the term with i = i 

t for vice versa) is 



1, m 



contributes, and from the sum ^^-^ only the term with / 
nonvanishing. Further, for a one-kink configuration t is always even. The result of the 
summations [^^^^^ Y^p] is again independent of i: 



kink 



2("-'H9^n3) ■ (S.ns) . 



(72) 



In order to carry out the sum over all possible one-kink configurations, we need the form of 
the function d'^Il^'^'^^ / dq'^dq\\q=Q for a general one-kink configuration of length n. If it starts 
with ri = (— ), we have 



2Y\kink 



dqpdq\' 



2-1 d^9o{k - q) 



dqt^dq'^ 



q=0 



(73) 



We add the configurations starting with ri = (— ) and ri = (+) and perform the sum 
Yl^=2teven ^^^^ possiblc ouc-kiuk Configurations. We find for the contribution of the 
one-kink configurations to the sum Yn r„ ^'i- O)- 



n{kink) 



2(n-l) 



5'o,+ (A;)2 1 - go,-ik)^^ 



-1 



^o,+(^)"^ -9o,-{k) 



-1 



d^9o{k - q) 



dq'^dq" 



g=0 



d^xdr {df,n^)-{dyn^) . 



(74) 
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The contribution of the mixed derivatives to Tr(GoS'^'^^)" from Eq. (|43|) is 



n-l 



dq^dq: 



qi,j=o 



ij ' 



(75) 



where T^j is given by Eq. (^2|). As we will show, nonvanishing contributions to the mixed 
derivatives originate from the strictly alternating configurations of the {rj} as well as from 
configurations where the alternation is broken up to two times. We will again study the case 
of the alternating configurations first. As can be seen from Eq. (^), Tf^ falls into two parts, 
the sum J2lm=i X]i=i Sm=j+i- '^^^ result of the summations [X]«m=i S'p] 

independent of the value that j takes and its contribution is given by % in Eq. (|69D. The 
terms occurring in the sum ELi Ei=i+i are Tj^il^^^.y ^(^=,r^±l) and Tj^il^^^z+iy "hose 
contributions are given in Eqs. ([66|), (]67|), and (pS]). 

We determine the number M of each of these three types of terms in the sum 
Yl\=i Sm=i+i every value of i and j, and denote the result of the summa- 
tions [{Y!i,ra=i + T!i=iY.L=i+i)T^p] by i-e. %j 



M(z=^±i) T(';^^^i) + M(i=„±2.+i) Tl^Lm±2z+iy We find 



7-,(z = odd, 3 = even) = -2^''-^\d^n^) ■ {d,n^) , 
even, j = odd) = -2("-3)(9^n3) ■ (S.ng) , 
even, j = even) = -2("~^)(9^n3) ■ (S^ng) , 

2 

7^,(z = odd, 3 = odd) = 2("-3)(9^n3) ■ (S^ng) - 2("-2)^(9^n,) ■ {d,n,) . (76) 



b=i 



In order to perform the summation 'Yllj=i i^j^ determine the second derivatives of the 
product of propagators from Eq. (|39|) for an alternating configuration starting with ri = (— ), 



dqidq"^ 

■^''odd ■^Je.ve.n 



qi,j=o 



J2 aoik)"^'' ^oAk)"^-' 

k 

dgoik-q) dgo-{k-q) 



dq^ dq\ 



even -^Jeven 

2rr-+ 



dqLdq 



^odd 3odd 



k 

goik)"^ go,-{k)^' 



dqt" 
n dgo{k - q) 



dq^" 

2 dga^k - q) 



q=o dq" 
dgoik - q) 

g=0 



g=0 



g=0 



^q^' 



dq"" 



dgoAk - q) 



g=0 



dq" 



q=0 



(77) 



We add the alternating configuration starting with ri = (+), and obtain from the summ- 
mation Xirjii i^j contribution of the alternating configurations to the sum Ylin r„ 
Eq. (0) 
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2) 



k 



dgo{k-q) dgo{k-q) 



^q^' 



q=0 



(7=0 



dgo-{k-q) dgo_{k-q) 



^q^' 



g=0 



dq" 



q=0. 



d^xdr q - 1) 5^(9,n,) ■ (d^n,) - - 1) {d,n,) ■ {d.n. 



b=l 



+ [J29o{k)^~'9o,-{k)'^- 

k 



n-idgoik-q) dgo^k-q) 



dqf^ 



q=0 



dq'^ 



q=0. 



X 



[Tl 
n(- - 1) {d^ns) ■ (duUs) 



Next we discuss the contribution of the nonahernating configurations to the mixed deriva- 



and T^iLm±2z) vanish. As in the 



tives. For nonalternating configurations the terms Tg^^^^^ 
case of the quadratic derivatives, the terms T^^^^j^^-^ and Tl^iLm±2z+i) nonzero if the 
contraction of / and m restores the alternation of the remaining sequence. For the mixed 
derivatives considered now, this is possible if the alternation of the configuration is broken 
not more than two times. We will first evaluate the case of the one-kink configurations. 

The first sum in V^^ (see Eq. (^)), has been discussed for one- kink configurations 

already, so that we just recall the result for T^"^^^^ Eq. ([7^). In order to evaluate the sum 
'Yl\=i'Yl,m=i+ii fi^^^ take i < j. The number of the site after which the kink occurs 
is again denoted by i. A contraction of the derivatives can restore the alternation of the 
remaining sequence only if I = 1 A m = t or / = lAm=(6+l). The sum ^""i^ ^-^j only 
yields nonzero contributions for every i with i < j if j = i (where only the term with 
/ = 1 A m = rrimax = ^ contributes to Yll=i X]m=i+i X]p)> ^^r every j with j > i if i = t 
(where only the terms with / = 1 A m = nimin = (^ + 1) contributes to X]l=i Sm=i+i J2p)- 
The summations [{J2l^m=i + SLi J2ln=i+i) J2'p] '^^^ be calculated to give 

t) = 2^^-'\d,ns) ■ {d,n,) , (79) 
t,j>t) = 2^^-^\d^ns)-{d,n,). (80) 



kink I 
ij \ 



At the present step of the evaluation, the contribution of the one-kink configurations starting 
with ri = (— ) to the mixed derivatives takes the form: 



g^n 
2, 



_2(n-3)^ 



n-2 E-1 



EE 



2Y[kink 



t=2 
ieven 



q=0 



n-l 

+ E 

j=t+i 



2Y\kink 



dq'^dq'' 



q=0. 



j d^xdr((9^n3) ■ {d^ris] 



^1) 



As to the second derivatives of the product of propagators from Eq. (|3 
are encountered: 



the following cases 
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kink 



9of do' 



2T\kink 



kink 



J2 9o{^y ^ 9o,-ik)hoAky 

k 

dgoik - q) dgo{k - q) 



dqi" 



9=0 



dq'' 



q=0 



dm 



dq'^ dq 



^go{k)2 ^go^ky^ ^goAkf' ' 



dgo{k-q) dgo-{k-q) 



dq'' 



q=0 



dq" 



9=0 



kink 



dqtidql^. 



/_^goiky^ ^ go,-ik)'^ go,+ {f^)^ ^ 

k 

dgoik -q) dgo^k-q) 



<li,j=0 



dqt" 



q=0 



dq" 



9=0 



^2) 



In the summations over i and j of Eq. ([81|), terms with [i even, j = t) or [i = i, j even) 
appear (| — 2) times, terms with {i odd, j = t) appear | times, and terms with {i = t, j odd) 
appear (f — |) times. Accounting for the case i > j and adding the one-kink configurations 
starting with ri = (+), we finally perform the summation ^ 
for the contribution of the one-kink configurations to the sum Yin 
mixed derivatives: 



n-2 

i=2, ieven 



and arrive at the result 



in Eq. (|75| ) for the 



n{kink) ' 



n 2("~2) 



go,+ {kyi 1 - go_{k)2 



--1 



|(n-4)^o,-(A;) go,+ ik) 



[go,-ik) - goAk)] 
dgoik - q) 



X 



+ ^g^{k) 



dq'' 

d[go,-{k - q)goAk - q)] 



dq" 



+ Yl Soik) 



1 - 



goAk) 



,n 



(--2) [goAky 



n 



2\9oAk)goAk)'' 



goAk) 

goAk)goAk)'^'^ 



-1 



X 



dgoik -q) d[goAk-q) ^goAk-q)] 



dq'' 



q=0 



dq" 



q=0 



q=0 
9=0- 

- goAk)' 



dx{d^,ny)-{dyny) . (83) 



What remains to be evaluated is the contribution of the two-kink configurations to the mixed 
derivatives Eq. ([75l). The first sum in F^^, Yl\m=i^ gives zero for configurations where the 
alternation is broken more than once. Thus we are left only with the sum Yl\=i Sm=i+i' 
where we first take i < j. The number of the sites after which the first and second kink 
occur are denoted by and I, respectively. Note that t, I are either both even or both odd. 
The sum YliiJ=i i^j only yields a nonzero contribution for i = i A j = I, where only the term 
with / = Imax = B V m = rrimax = I contributes to J2m=i+i Yl'p] ■ Thus the summations 



m=j+l 



E'p] yield 
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1-3 



two— kink 



(84) 



Upon the determination of the second derivatives of the function 11 for a general two-kink 
configuration, it turns out that after adding the configurations starting with ri = (— ) and 
ri = (+), the form of d'^Il^'^"'^™^ / dq^^^dq^^^^ depends only on the distance between the 
kinks. Introducing c) = [ — 6, we have 

^2 j~jf«)o— femfc 



qi,j=o 



dgojk - q) dgojk - q) 
dq^^ q=o dq'^ 



g=0 



15) 



Now we carry out the sum over all possible two-kink configurations by summing over all 
possible distances D {d runs from 2 to {n—2) over all even values) and weighting each distance 
by a factor {n — d — 1), which is the number of different two- kink configurations with the 
same inter-kink distance. We find the following contribution of the two-kink configurations 
to the sum Er-i...r„: 

n , goAk)^'^ - 9oAky~^ 



v-ij 

n(two—kink) 



2(n-2) 



dgo{k-q) dgo{k-q) 



dqt" 



q=0 



dq" 



q=0 



goAk)-' - goAk)-' 



n-2) d'^xdr{d^n3)-{d^n-i). 



^6) 



In order to obtain a closed expression for the fermion determinant, we still have to re-sum 
the power series in n from Eq. (^) to infinite order, namely 



1 



V-Tr(GoS 
^-^ n 



(0)^ 



n=l 



E X^^ -u Y^^ -u Y^J -u Y^J -u Y^J 
\^n{alt) "T ^n(kink) ^n{alt) ^n{kink) ^n(two-kink) 

n=l 



17) 



where the terms on the right-hand side were defined in Eqs. (|7TD, (0), (0), (P5|), and (^). 
The summation is achieved in terms of the polylogarithmic functions Lio, Lii, and Li2, 
which have the power series expansions Lio (2) = -2^"; Lii(z) = J2'^=i ^ — ^^T^l' 

Li2(-2) = Yl'^=i From the integral representation of Li2{z) it follows, for example, that 
^Li2(-2) = —liKlzfl^ This property will be needed when the derivatives of these functic 
as shown in Eqs. (El 



tions, 



pop below, are to be performed for a later numerical evaluation. In 



addition, we use the relation 



.„=i nz - 

the contribution of the second-order derivatives to the helical part, Eq. 
determinant: 



. After some more algebra, we obtain for 

of the fermion 



n—2 

neven 



2 

d^a;dr ^ d^Ua - d^Ua + x 

a=l 
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where we have introduced the abbreviations 

1 ^2 



X 



1 92 



q=0 



8a2 dqi^dq'^ . 
1 



g=0 



+ 



4a2 



^9) 



(90) 



Here ^f{q) was defined in Eq. (p?]), and the exphcit expressions for the functions $f (g), 
$f (g), ($^)'''', ($6 )^'' hsted in the Appendix. As can be seen from Eq. 

and Eqs. (|A1|) - (|A5D , the parameters of the initial model enter in a nonperturbative way into 
the functions $, which describe the infiuence of the doping. 



B. Ferromagnetic contribution to the fermion determinant 

We are left with the task of evaluating the part of the fermionic determinant which 
contains the ferromagnetic fiuctuations L, namely the second trace on the right-hand side 
of Eq. (|3TD. Again we expand up to second order in powers of a. Whereas in the spiral 
contribution TTln{GQ^ — S^^-*) every order in the expansion of the logarithm had to be kept, 
the ferromagnetic contribution directly generates an expansion in powers of a: 

Trln[l - aG'oS(^) - a^GoS^^)] = Tr [-aGoJ^^'^ - a^GoS^^) - ^GoS^^^GoS^^)]. (91) 

In order to evaluate this trace, we must determine the inverse of G^^ = Gq^ — T,^^\ While 
Gq^ is diagonal in spin and momentum indices (see Eq. (0)), S'^^) splits into two parts 
depending on the wave vectors {k — k' — kg) and {k — k' + ks), respectively (see Eq. (^71) ). 
Since we consider long-wavelength spin fields, {k — k' ± kg) must be small. Thus, S^^^ has 
matrix elements only on four side diagonals between k and k±ks. Hence Gq ^ has the form 
of a band matrix, where the number of matrix elements is determined by the number of k 
points and every k point splits into a 2 x 2 spin space. Obviously, it is impossible to invert 
Gq ^ exactly. However, Go can be found perturbatively in powers of a by finding a matrix Go 
which, when multiplied on Gg^, gives a diagonal matrix V plus off-diagonal terms of 0{a): 

Go ' Go = V + Oia) . (92) 

This works because the square of the matrix S*^*^) is proportional to the unit matrix 
(S^°^)L,a'fc' = a^Saa'hk'. Based on Eq. (P, we define (si°))^,^^,,, = laV^{k - k' ± h), so 
that = SL°^ + Sf . Now we try as an ansatz for Gq: 

Go = Pi + P2 + , (93) 
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multiply it onto Gq^ from the right and then choose Pi and P2 so that Eq. (p^) can be 
fulfilled. An appropriate choice of Pi, P2 should make terms which are manifestly non- 
diagonal in /c-space either zero or small and of 0{a). Moreover, terms which were already 
diagonal in fc-space should not lose that property through multiplication by Pi, 2- This is 
accomplished by choosing 



Pi = Sf Gq-L /g' 

p, = sf G^X 1 9' 



(94) 



where Gq _ and Gq _,_ are exactly diagonal in spin and momentum space. Their /c-dependence 
is chosen to be 



(<^0,±)afc,a'fc' - fl'0,±(^) Saa'Skk' 



(95) 



Then from the multiplication 



(G, 



-1 v(o) v(o)^ 



EV^O(^i + ^2 + + ) = D-F + K_ + K^ 



(96) 



we obtain a purely diagonal contribution 



(97) 



a contribution being approximately diagonal in momentum space but off-diagonal in spin 

space 

iF)^k,a'k' = \go\k)al^,nl{k-k')Q_{k'), (98) 
and terms which are off-diagonal in spin and momentum space, but which are of 0{a): 

K± = Go - si^^Go-^ 



Mk) - go{k' T ks)K^,n''_{k - k' ± k^) . 



(99) 



We have used the abbreviation: 



Q±{k) = g,\{k) ± g,X{k) . 



(100) 



In order to diagonalize the remaining problem, we multiply Eq. (|96D by a matrix (D + F) 
from the right, where 



(D)afc,a'fc' = Q-{k) ^ g^ ^{k) g+{k) - g"^ 



^aa'^kk' , 



{F)ak,a'k' = - gZ\k) n^(/c - k') g,\k') 



(101) 
(102) 
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This leads to the desired expression Eq. (|92D. The matrix Go which we have been looking 
for is thus given by 



Go 



^(0) ^(0) ^^ljg2 ^ g-2j.i0) ^(0) ^^ljg2 ^ j.,^> ^ J. 



(0) , v(0)' 



D + F 



(103) 



(104) 



(105) 



(106) 



(107) 



and for the diagonal matrix V in Eq. (0) we obtain 

{V)ak,a'k' = DD - FF = 

= Q-J{k) [g' - g' g^\k) g^{k) + gf{k) g,^_{k) g^Xi^) 
In order to simplify the notation, we define 

d{k) =g'-g' go\k) g^{k) + gf{k) g,l{k) g^Xi^) ■ 
The small off-diagonal terms are 

aR = (DF - FD) + {K_ + K+){I) + F) . 
From Eq. ([9^), the inverse of Gq ^ to first order in a is then found to be 

Go = Gol?"^ - aGoV-^RV~^ . 



This has to be inserted into Eq. (0), together with the expressions Eqs. (|2^), ( p9| ) for 
the first and second orders of the fermionic self-energy. After inserting Go from Eq. ( |103| ), 
Tr[aGoS*^^^] then becomes a sum of 240 terms, while Tr[a^GoS'^^)] contains 112 terms and 
Tr[^|-GoS(i)GoS(i)] becomes a sum of 1600 terms. (We classify different terms by differ- 
ent combinations of the order parameter fields. Every individual term still comprises the 
momentum- and spin-space summations arising from the matrix multiplications). These 
numbers grow considerably when the strings of Pauli matrices contained in these traces are 
evaluated, e.g. a string of six Pauli matrices leads to a sum of 15 different permutations of 
contractions. While strings of even numbers of Pauli matrices are carried out using the trace 
reduction formula Eq. (p6D, strings of odd numbers of Pauli matrices are evaluated by replac- 
ing two Pauli matrices according to the identity J2a2 '^a\a2'^ala3 = ^"^"^^aiaa + i£''^''^''''^o'^{a3- 
Obviously, a string of an odd number of Pauli matrices leads to inner products and vector 
products of the involved vectors. The traces are evaluated with an algebraic program^ which 
we developed using Mathematica It reduces the trace over the Pauli matrices, performs 
the contractions of the vectors, implements the constraint Eq. ([19D , and takes the trace in 



The program TraceEval is available via anonymous ftp from ftp.physik.uni-wuerzburg.de as 
/pub/dissertation/klee/TraceEval.m. 
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momentum space by keeping only terms which are diagonal in the momentum indices. The 
result for the ferromagnetic contribution to the fermion determinant is 



ferro 



J d^xdr jxi + X2 



[L-n.y + iL-n^) 



+ X3 (L-Usf 



a 



[L-ni){ni-d^n3) + (L-n2)(n2-(9^n3) 



2X4 r n ^ , Xg 



Xe 
a 



L-df^ns H {L-n-i){ni-d^n2) 



[L-ni){n2-d^n3) - {L-n2){ni-d^ns) 



+ L-{ni X d^ni + n2 x d^n2j + ( ns x d^n^ 



:i08) 



The quantities Xi • • • Xs; X4 • • • Xs generalized susceptibilities of the fermions in the pres- 
ence of the spin fields. The explicit expressions are listed in the Appendix. The functions x 
consist of summations over e„ and k of polynomials of the inverse free fermionic propagator 
and its derivative. Shifting the summation index k —k and using the symmetry of 
the free propagator, Eq. (0), it can be shown that 



X4 
X5 





Xe 



A* 
X7 



Xs 







for 
for 



/i 



0, 
1,2, 



(109) 



i.e. Xa is non-zero only for the spatial components /i = 1, 2, while Xs ; Xe ; X75 Xs non- 
zero only for the time-component /i = 0. Note that the susceptibilities are proportional to 
d~^{k) or d~'^{k), where d{k) was defined in Eq. ( |105| ). The zeroes of d{k) determine the 
dispersion relations of the holes in the helical spin background. Since d{k) is fourth order in 
the free inverse propagators g^^, which contain contributions from two bands, d{k) is eighth 
order in ie„, and eight hole bands will result. This additional splitting in comparison to the 



antiferromagnetic case, in which four hole bands were obtained [^, results from the fact 
that the momentum transfer on the fermions interacting with the spin field is now Q ± 6Q 
instead of Q as in the antiferromagnetic case. 



C. Contribution of the Heisenberg part 

Finally we have to calculate the continuum limit of the pure spin part of the action, 
Eq. (^. The study of the continuum limit of this action, to which our model simplifies in the 
absence of doping, has in itself attracted much interest. In order to investigate the effect of 
frustration on the stability of the ordered state against quantum fluctuations, field-theoretic 
mappings have been performed for the Heisenberg antiferromagnet on a triangular lattice 
1^-^, for frustrated Heisenberg antiferromagnets on a chain and on (i- dimensional lattices 
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36| , |40| . The frustrating effect of either the lattice geometry or the competition between 
nearest-neighbour and higher-order interactions leads to noncollinear classical ground states, 
so that the order parameter is an element of SO (3). Thus our gradient expansion of the 
pure spin part, carried out in position space, is similar to the above-mentioned treatments. 
However, in these works the spins in the ground state have a periodicity, so that one can deal 
with a finite number of sublattices. This does not apply in the case of the incommensurate 
spiral. 

It is useful to extract the antiferromagnetic modulation from n by writing Upg = 
(^—l^ip+g) nipg, where we define 

"^pg = 'Til cos{SQ-rpg) — n2 sm{SQ-rpg) (HO) 



First, we consider the Berry phase. We inject Eq. (|2^ for Spg and expand A{Spq/S) and 
drSpg to first order in a. We exploit the gauge freedom that we have in the definition of 
A to choose it an even or an odd function of its argument, use the constraint satisfied by 
the monopole potential, Eq. (^), and the fact that total time derivatives drop out upon 
r-integration. Discarding terms which are of third or higher order in derivatives of m, we 
obtain: 

-^Berry = / dr i5 V [(- 1) A(mpg) ■ S^TTipg + aL- {rupg X drinripg)] . (Ill) 

From the term involving the vector potential in Eq. ( |1 1 1]) , we first consider one chain of 
lattice sites by holding the index q fixed: 

''PQ 



SA{q)= [ drz5V(-l)^A(mp,)-a.m, 
Jo p 

= j driS" ^ [-4(mpg) ■ drinripg - A{mp+i^q) ■ drfUp+i^g . (112) 



Expanding our smooth order parameter fields according to ni^2(p +1,?) = i^i,2ip,Q) + 
cidxni^2{p,(l), we find: 

TTip+ig = mpgCos{6Qxa) — rhpgSm{6QxCi) 

+a (A^.mpq) cos((5(5a-a) - a {Aj^rhpg) sin{6Qr,a) , (113) 

where we defined rhpg = rii sm{6Q-rpg) + n2 cos{6Q-rpg). Note that in contrast to collinear 
spin configurations, the zeroth order contribution to rUp^i g is not rripg, but a vector which 
we denote by Tn^^^^^ = mpgCos{6Qxa) — rhpgSm^SQxCi)- With expressions like A^rripq we 
mean: 

A^mpg = (d^ni) cos{5Q-rpg) - {d^n2) sm{6Q-rpg) . (114) 
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After some algebra one obtains 



(115) 

Now, the terms dTA{mpg) ■ drrripg and /^^ drA(mp_,_j^ g) ■ 9^mp_,_j^ g describe the area of 
the caps bounded by the paths of rripq and mp_^-^^ g, respectively. Since nipg and rrip^^^ g are 
connected by a space-time independent rotation about the axis by the angle (SQ^a), the 
solid angles spanned by their paths are equal. Thus the corresponding terms cancel when 
(|115|) is inserted into Eq. ( [L12| ). After some tedious algebra, we get for the remaining term: 

mj_^i g ■ {A^ml^^ g X A^mJ+i J = -nr(9^ni x 9^ni) cos{6Q-rpg + 6Q^) 

+n2-{drn2 X dxn2) sm{6Q-rpg + 6Qx) (116) 

Employing the relations e'^^'^n'^drn^^Oxn'^ = ^defdr'n'^ d^nj, this result can be shown to agree 
with the one obtained for a spiral configuration in a frustrated Heisenberg chain within a 
different method Due to dimensional reasons, we must have: SQ = Sip /a. We assume 
that in the continuum limit the angle 6cp between neighbouring spins remains constant. Thus 
in the limit a —>■ the terms containing cos{6ip-x/a), sm{6ip-x/a) oscillate very rapidly and 
drop out. Then the total Berry phase is given by the second term in Eq. ( |1 1 1|) , from which we 
keep only non-oscillatory contributions. Taking the continuum limit lima^o Ylt 
leads to 



pq 



Jd^X 



We now turn to the real term in Eq. 
1 



iS 

dr / d^x — L ■ {rii x drUi + n2 x (9^-722) 



2a 



(117) 



). We attribute to each site the quantity 



M=±l S'=±l ^ 

and then perform the sum over all sites. As before, we insert Eq. ( P^ ) for S and expand the 
fields at site {p+6, q) in terms of the corresponding fields at site (p, q) using Eq. ( |113|) . Again, 
we eliminate short-range oscillatory contributions. After some lengthy but straightforward 
algebra we obtain 

/3 . .2 



5, 



real 



- / dr 



d^x Jh S"^ < a ^ sm.{5Q^a) rii ■ d^n2 



+ ^-cos{SQ^a) {df^niY + {d^n2] 



1 2 ^2 

2 + - ^ cos{5Q^a) - ^ ^ cos^(5Q^a) 



y^cos(^Q^a) 



(L-mf + (L-n2) 



(119) 



Note that here the index /i runs only over the space components, namely /i = 1,2. 
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V. DISCUSSION OF THE CONTINUUM THEORY 



In order to express our continuum action in terms of the SO (3) order parameter Q = 
(ni 7x2/13), the ferromagnetic fluctuations L have to be integrated out. Since our resulting 
action is linear and quadratic in L, this can be done by extremizing the action with respect 
to L. Inserting Eqs. (|108|) , ( |117] ), and ( |119| ) into the saddle-point equation -^{S{evm + 



5, 



Berry 



real I 



0, we obtain 

[Xi + X2 {nln\ + nln\) + Xs = (120) 

where v denotes the sum of all terms which were multiplied by L linearly in the action 



1 

2^ 



i=l 

2 2 



:i2ii 



i=l 



and we have introduced the following abbreviations to incorporate fermionic and Heisenberg 
contributions: 



Xi = Xi- JhS^ 



2 + ELi cos{5Q^a) 



X2 = X2 + JhS^ 2 + -ELi cos{5Q^a) - -ELi cos^(^Q 



-At=0 
X7 



/t=o ^"S" 

x^ -Y 



(122) 



We solve the saddle-point equation for L by multiplying both sides of Eq. (|120|) by the 
inverse of the matrix [xi (^"^ + X2 {nin\ + ^2^2) + Xs ^3^3] ^^id obtain 



-<[C2"n?a„n^-C3°n^M 
-nl[Qnld^nl + Qn\drnl] 



(123) 



The newly introduced abbreviations Ci, C2 ; C3 ? ^-^id C4 are defined in the Appendix. Recall 
that the upper index on the generalized fermionic susceptibilities denotes their frequency 
component, which accompanies the time component of the field derivatives, while the index 
a = 1,2 runs over the wave- vector components of the susceptibilities, which accompany 
the spatial components of the field derivatives. In the above solution for Eq. ( |123| ), we 
have applied the results of Eq. (|109|) to show explicitly how L is composed of first-order 
space and time derivatives of Q = {ni n2 n^). Our expression for L must now be reinserted 
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into Eqs. ( |108| ), ( |117|) , ( |119| ) for 5'forro, 5'Berry, and S'reai- In order to simplify the result, 
the following identities, derived from the constraint on the order parameter, Eq. ([T9|) , are 
needed: 



(ni ■ 9^.7x3) (n2- (9^ ng) - {ni-d^n3){n2-df,n3) = ■ (S^ng x d^n^) 
d^ris ■ (ni X d^,ni + x S^na) = ng ■ (9^ rig x d^n^) 

The total -L-dependent part of the action can then be cast in the form 

Sl = -X^" j d^xdr ((9^ni ■ d^Ui + S^ns ■ S^ng) 



(124) 



- (r^ - X 



X 



/d^xdr„,.(a,n,x9„n3) 



where the prefactors are given by 

1 



2a2 



[Xl + X3 



-1 



^ [Xi + X2] ^ [X4X4 + Xlxl + X7X7 + XsXs 



(125) 

(126) 

(127) 
(128) 



From Eqs. (|60|), (|88| ) - ( pOl) for the spiral contribution to the fermion determinant, Eq. ( 119 ) 
for the contribution of the Heisenberg part, and Eq. ( |125| ) for the contribution of the in- 
dependent part of the action, we find the total continuum theory to be given by an SO (3) 
quantum nonlinear a model, a term linear in the derivatives and a geometric term which is 
third order in the fields and contains first-order derivatives with respect to time and space. 



The final form of our quantum nonlinear a model action is 

3 

'5'qnlo-m 



^ a=l 



X] 



(129) 



(130) 



where Pi" = p^ 7^ "p^ . The coupling constants of the model are given as a function of the 
microscopic parameters and the generalized fermionic susceptibilities: 



X 



cos(5g^a) 6^"- 



(131) 
(132) 
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Equivalently to Eq. ( |130|) , the model can be written in SO (3) matrix form: 



(133) 



iiv 
Pi 





av 
Pi 



\ 



; 

taking into account the constraint Eq. 



(134) 



there are 



'qnlcrm 

where the coefficient matrix is given by 

/ 

pt^u = 

V 

Among the nine fields Qab{.x) = nl{x) 
three independent fiuctuating fields. Hence, our continuum action describes three massless 
modes, or spin waves, which result from the complete breaking of the 0(3) rotation group by 
the ground state. Each spin wave excitation corresponds to infinitesimal rotations about one 
of the n^'s. One Goldstone mode describes variations of the spin orientation within the plane 
defined by the classical ground state, and thus corresponds to infinitesimal rotations about 
the axis n^, while the other two Goldstone modes describe out-of-plane fluctuations, and 
are related to rotations about rii and n2. Eqs. ( |130| ) or ( p,33| ) describe the linear spectrum 
and the interactions of these spin-wave excitations. 

The coefficient matrix P in Eq. (|134|) is represented in internal space. Every matrix 
element of P^*^, however, is still a matrix in the space-time indices fi and u. We will now 
consider the symmetry properties of the matrices Pi'^ and p'^'^ . As regards Pi'^, it can be seen 
from Eq. ( |131|) that it is symmetric in fi, v when looking at the definitions of i Eq- (H), 
and 1 Eq. ( p.2(j| ), and taking into account Eq. ( |109D , which implies that x!^^ is nonzero 
only for /i = z/ = 0. Inspection of the terms appearing in Eq. ( |132| ) for p^ (see Eqs. (pO|), 
(|126|) and ( [L27| )), shows that p^ contains also antisymmetric contributions. However, since 
we sum over /i and v in Eq. ( |13(J| ) or ( |133| ), and the term {dfj_na{x) ■ di,na{x)) is symmetric 
in fi, u, the antisymmetric contributions cancel and we are left with a symmetric coefficient 
matrix. 

Let us now consider the matrix elements P^'^ which mix frequency and wave-vector 
indices (or, in the corresponding field derivatives, space and time indices). For the spiral 
contributions to the generalized fermionic susceptibilities, it can be shown by performing 
the derivatives and using the symmetry of the free propagator, Eq. ([T^), that for a = 1, 2 

'^l,2,3(?) 



QqoQqc 



q=0 

and = X°° 



QqoQqo 

(•^4^,5,6) 



9=0 



(135) 



(136) 



As a consequence, x"" = ^ind x°°' = x^° vanish, as can be seen from Eqs. (^) and (pO|). 
For the ferromagnetic contributions, we know that x!^^ ^ only for fi = u = and that x^'' 
does also not mix frequency and wave-vector indices (see Eqs. ( |109| ) and (|l27D), so that 
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Pi 

P3 



Pi 
P3 



0. 



(137) 



Therefore, the matrices p'^'^ and p'^'^ are of block-diagonal form (for clarity, we will use r 
instead of for the frequency component of the coefficients, and a = x,y instead of a = 1, 2 
for their wave- vector components). 



/ p7 \ 



Pa 



pT p7 



1,3. 



:i38) 



V pT Pf I 
The action, Eq. (|130|) , can then be written as 



qnltrm 



d xdrTr 



a,f3=x,y 



(139) 



This form, in which the coefficient of the time derivative term is a tensor in spin space 
and the coefficient of the space derivative term is a tensor both in spin space and in real 
space, was suggested as the most general case for systems with noncoUinear spin orientation 
within a hydrodynamical theory |^1]]. For some special cases of the direction of the spiral 



wave vector, the matrices Pa become diagonal. If, as suggested by experiments |jT9| on 
cuprate superconductors and by theoretical studies |^ of models related to the spin-fermion 
model, the spiral wave vector lies along the (1,0) or (0,1) directions of the lattice, i.e. kg = 
(^ ± 6Qx, ^) or kg = (^, ± 6Qy), then, using the symmetry relation Eq. ([T5|) it can be 
shown that 



dq^dqi 



g=0 



1^4,5,6^ 



X 



xy 



0, 



0. 



(140) 



Thus, Eqs. (|T3T]) and (|T3|) yield 

(p7 \ 
Pa= . (141) 

V pyy J 

In the general case of Eq. (|138|) , the symmetric (2 x 2) matrix in x and y can be diagonalized, 
which amounts to changing to a new basis which is a linear combination of the lattice basis 
vectors. 

In order to make the properties of the spin waves resulting from our theory more trans- 
parent, we write the action in another equivalent form 
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Vi^m = -J d'xdr Tr [PZ, iQ^d.Qf] (142) 

(where we assume the coefficient matrix to be diagonahzed) . Now we represent the three 
degrees of freedom contained in the antisymmetric spin-space matrix {Q'^d^Q) according 
to ((5^^5^jQ)i2 = (^^^^/»^)i3 — -^/i; ^'^d ((5^^5^<5)23 = ^^l■> which leads to the action 

Vi.m = - j d^dr [2pYAl + {vV + V7){BI + Cl) 

+ E [2pr^' + +pr)(i?« + c^')]} (143) 

a=x,y 

The dispersion relations for the three spin-waves can now be read off: 

Pi Pi 

Thus, the spin-waves have different velocities which may depend on the direction of propa- 
gation, as was suggested by the hydrodynamical theory [^. It is the fact that the theory is 
not Lorentz invariant, i.e. that the coupling-constant matrix P^^'^ for the time components is 
not proportional to that for the space components, which allows the three spin-waves to have 
different velocities. In our case, the two out-of-plane modes have the same spin-wave veloc- 
ity, which differs from that of the in-plane mode. In the case in which the coupling-constant 
matrix is isotropic, i.e. P^^ = P^^, the spin-wave velocities become 

Ca = VpT/p7 (145) 
in the plane spanned by the classical ground state, and 

cb,c = VipT+pr)/ipr+pv) (146) 

out of this plane. 

Now we briefly consider the internal-space symmetries of the action which we have ob- 
tained in Eq. (|130|) or ( |133| ). This action is invariant under global left 0(3) rotations 



Q —>■ UQ, where U eO(3). This symmetry corresponds to the usual invariance under rota- 
tions of the basis vectors in spin space: — > U"''^ n\. In addition, the action is invariant 
under global right transformations Q QV if [P, t^] = 0. The dimension of the group 
to which V belongs depends on the values of the three coupling constants contained in P. 
In our case, in which two of the coupling constants are equal, VeO{2). This right trans- 
formation on Q corresponds to a mixing of the basis vectors: ra^ '^b^t^bc- Since 
Pi = P2 7^ P3; the mixing occurs between the basis vectors ni and n2. The right 0(2) 
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invariance reflects the screw-axis-like symmetry of the spiral state, which is invariant under 
a rotation about the axis by an angle enclosed between a pair of spins, followed by a 
lattice translation along the direction which connects the two spins. Thus, our action for 
the quantum nonlinear a model possesses an 0(3)®0(2) symmetry as in previous work for 
frustrated quantum spin systems jOHTT . 



In addition to the SO (3) quantum nonlinear a model, we find in our microscopic deriva- 
tion the contribution (see Eqs. (|60D, ( |119| )) 

S\in = y^A^ TT— "^^(5) — sm{6Qaa) / d^xdr rii ■ daUi , (147) 

a=x,y 

which is linear in the derivatives. The action 5*1111 is of course invariant under rotations of 
the basis vectors, n'^ U"'^ n^. Under the transformations Sfe'^ft He, however, S'un is 

invariant only if det = +1. Transformations with det V = —1 change our right-handed set 
of basis vectors into a left-handed one. This means that S\\n is not invariant under a change 
of helicity of the spiral. 

Since the linear term is not positive definite, the weight of some field configurations in 
the path integral will tend to infinity, and hence this term leads to instabilities. To recover 
a stable ground state, we must ensure that the action is at a minimum. Thus, we impose 
the condition that the prefactors of the linear contributions must vanish: 

1 d 



a 



<l>r(g) = Jh^' sin(5Q,a) . (148) 



g=0 



This condition yields two equations (one for a = x and one for a = y) for the two spiral pitch 
parameters 5Qx and 5Qy, whose values are then determined as a function of the microscopic 
parameters and the doping concentration. Thus, the wave vector of the spiral background is 
self-consistently determined from the stability argument. In order to check the consistency 
of the stability condition, Eq. ( |148| ) , we consider again the special cases of the spiral wave 
vector mentioned above. For a spiral wave vector in the (1,0) direction, the condition in 
Eq. (|148|) is trivially fulfilled for a = y because both sides vanish (the symmetry relation, 
Eq. (p!5D, implies that -^^i{q)\q=Q = for this choice of ks). If the spiral wave vector lies 
along the (0,1) direction, the terms on both sides of Eq. ( |148| ) vanish for a = x. All other 



choices for kg lead to nonvanishing terms on both sides of Eq. (|148|) for a = x,y. 

Finally, our microscopic derivation gives rise to the following term which is third order 
in the fields and second order in the derivatives: 



(149) 



From the definition x'"' = ^ [Xi + X2] [Xe + X7 + X8]xl (see Eq. (|T28D), and using 
Eq. (|109|) , it is seen that Sgeom is nonzero only for yU = r and u = x,y. For 6Q = 0, 
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can be shown to be zero, so that S'geom vanishes. Performing the first-order derivatives 
in Eqs. ([X9|), ( |A11|) - ([A13|) and looking at Eq. (|122|) , one finds that is purely imaginary. 



As is well-known, a term having the same structure as S'geom appears in the continuum 
theory of one- dimensional quantum Heisenberg antiferromagnets with an order parameter 
in 5*2. The Pontryagin index 

Q=— dxdr n ■ [dxTi X dru) (150) 
47r J 

describes the winding number of the mapping 712 (6*2) = Z. The long-wavelength action 
of the spin-S* antiferromagnetic Heisenberg chain contains the term = i27[SQ, which 



has led Haldane ||3^ to conjecture that integer spin chains posses a gap in the excita- 
tion spectrum in contrast to half-integer ones. In a two-dimensional square lattice, how- 
ever, the topological term is cancelled, since the summation of all rows of spins yields 



lima_»o ^27rS' Q(g) = / dymSdyQiy) = A similar situation arises in the 

triangular lattice, where the corresponding topological term does not contribute to the dy- 
|3[. Using the symmetry properties of the coefficient x'^'^, we can write our geometric 



namics 
term as 



-S'geom = Y.a X'"^ j "^3 " {dan^ X drU-i) 

= 47r x"" j dy Q{y) + 47r x"^ j dx Q{x) . (151) 

While Q is an integer, the coefficients in front of the topological terms 47rx^" are not nec- 
essarily integer multiples of tt as in the undoped models, but are in general functions of the 
microscopic parameters and the doping. Therefore, the obtained field theory contains two 



^-vacuum terms (one for each spatial direction) with continuously varying parameters, 
in contrast to the until now known field theories obtained for quantum spin systems in two 
dimensions. 

It should be noted that S'geom does not correspond to the topological invariant classifying 
the mapping 773 (50 (3)) = Z, which is third order in derivatives and therefore is not contained 



in our expansion. It is given by 



Q = l/(247r2) j d=^a;dre^^'^Tr [{Q'' dxQ){Q^ d^Q){Q^ d,Q)] , (152) 
If we write S'geom in terms of the S0(3) matrix Q, it becomes 

^geom = -r j d'xdrTr [C{Q^ d^Q){Q'^ d,Q)] , (153) 



where C12 = —C21 = 1 and Cij = otherwise, and we have introduced f^" = \{x^" ~ X 
which is completely antisymmetric in the momentum-space indices /i and v. A study up to 
0{a?') is left for the future. 
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Since the variation of 5'geom under an infinitesimal variation of the field vanishes, its value 
can only change in discrete steps and thus it is quantized. Field configurations falling in 
different homotopy sectors cannot be continuously deformed into one another. Non-trivial 
field configurations correspond in this case to instanton-like configurations for the plane of 
the spiral. Since until now the few studies carried out for the ^-term considered only the 
one-dimensional case such that an extension to (2+1) dimensions is not straightforward, 
we restrict ourselves for the following discussion to the sector in which S'geom = 0, such that 
the results obtained previously in renormalization group analysis can be applied. 

Renormalization group analyses of a (2 + l)-dimensional S0(3) quantum nonlinear a 
model in the form which we have derived in Eq. ( |142| ) have been performed by Azaria et al. 
and by Apel et al. 0. Although previous microscopic derivations of the SO (3) quantum 
nonlinear a model from frustrated Heisenberg models resulted in an action with pY = pY 
and pY = pf = the renormalization group calculations were carried out for the 

model which we obtained in Eq. (|142|) , because the aforementioned conditions are not stable 
under renormalization. Azaria et al. considered the case in which the coupling constant 
matrix is isotropic in lattice space, i.e. P^^ = P^^, so that the action contains four coupling 
constants, which correspond to pY , pY , pY i pT case. Instead of using these four 

coupling constants, it is possible to work with two independent spin stiffnesses, Pa oc pY and 
Pb oc {pY +pY)y two uniform spin susceptibilities, Xa oc pY and Xb oc {pY +pV)- (The 
spin- wave velocities are then given by = \/ Pa/Xa and Cg = a/ Pb/Xb-) In contrast to the 
case of the 0(3) quantum nonlinear a model describing the long- wavelength properties of 
collinear spin configurations, here the spin-wave velocities Ca and Cb already renormalize at 
one-loop order. In fact, the spin- wave velocity in the 0(3) quantum nonlinear a model at 
zero temperature is not renormalized at one-loop order [|T| as a consequence of the Lorentz 
invariance of the theory. Such an invariance is absent in the present theory, as already 
pointed out before. 

The one-loop renormalization group equations for the parameters c^, Cs, p^, p^, [3 admit 
a nontrivial fixed point for T = 0, c* = c* and p* = p*. At this point, the theory 
is 0(3) (8)0 (3) ~ 0(4) symmetric and Lorentz invariant. As in the case of the classical 
0(3)®0(2) nonlinear a model in (2 + e) dimensions |p!0| , pn^ (which contains two couplings 
in the isotropic case), the symmetry is dynamically enlarged at the fixed point. Connected 
to this fixed point in the five-dimensional parameter space is a critical surface separating an 
ordered phase, in our case a spiral one, from a phase disordered by quantum fluctuations. 
The scaling properties of this quantum transition have been described by Chubukov et al. 



12[| . The present work provides a way to determine quantitatively how a variation of doping 



will either drive the system from the helically ordered to the gapped spin-liquid phase, or 
in case the spin-wave velocities vanish, an instability to another state occurs. 
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VI. CONCLUSION 



In this work, the long wavelength, low-energy sector of a spiral spin configuration in 
the microscopic spin-fermion model has been mapped onto an effective field theory. This 
mapping has been accomplished by a systematic gradient expansion of the spin action which 
was obtained by exactly integrating out the fermions from a coherent-state path-integral 
representation. The magnetically ordered spiral state completely breaks the 0(3) rotation 
symmetry in spin space. This leads to an SO (3) order parameter. We have shown how the 
constituent fields of the order parameter are obtained from the physical spins, such that 
the ambiguity in describing the physical spin on 5*2 with elements in SO (3) is avoided. The 
low-lying modes fluctuating around the classical spiral ground state were parametrized by a 
decomposition of the spin fields into a helical and a uniform component. 

The gradient expansion of the fermionic determinant leads to an infinite series that can 
be summed to all orders of the microscopic coupling constants by a combinatorial method 
which exploits the constraint on the order parameter. We have carried out the gradient 
expansion up to second order, such that the relevant terms in the same order as those 
generally proposed in phenomenological approaches are contained. Such an expansion leads 
to an effective action containing first- and second-order space-time derivatives of the order 
parameter. 

The first-order terms yield through a stability condition two equations which determine 
the two spiral pitch parameters as a function of the microscopic parameters and the doping. 
Thus, our theory is able to determine the wave vector of the spiral spin-background in a 
self-contained way. 

The second-order terms, on the one hand, lead to an 0(3)(8)0(2) symmetric quantum 
nonlinear a model describing the long-wavelength spectrum of three spin waves. The co- 
efficient matrix of the a model, which is responsible for the Lorentz noninvariance of our 
theory, determines the spin-wave velocities and spin-wave stiffnesses as a function of the 
microscopic parameters and the doping. It should be emphasized that the continuum ap- 
proximation only refers to the spin dynamics, so that our method can deal with arbitrary 
dispersion relations for the fermions. Since the gradient expansion is nonperturbativc in the 
coupling constants of the microscopic Hamiltonian, its results are also relevant for related 
models, like the Kondo lattice model (J// = 0) and the t-J model (J/^ oo). 

On the other hand, our continuum theory yields geometric terms of the same form as the 
one obtained for the one dimensional antiferromagnetic Heisenberg model in the continuum 
limit. In our case however, the corresponding coefficient varies continuously as a function 
of doping. Whether these terms may close the gap in the quantum disordered phase for a 
given value of the coefficient (like at ^ = tt in nonlinear cr-models with a Q term in (1+1) 
dimensions) , is a question left for further studies. 
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Our results show that by deahng with a given but rather general model for doped an- 
tiferromagnets, new terms appear that are in general not contained in phenomenological 
approaches, where only symmetry arguments are used. Furthermore, the influence of the 
doped charge carriers on the critical behaviour of the incommensurate spin system is ob- 
tained explicitly. The doping dependence is contained in the parameters of our SO (3) 
nonlinear a model through generalized fermionic susceptibilities, which in addition yield the 
new dispersion relations of the holes in the presence of spiral spin flelds. The numerical 
evaluation of the generalized fermionic susceptibihties is presently being carried out. 

We would hke to thank A. Angelucci for useful discussions. S.K. gratefully acknowledges 
support by the Studienstiftung des deutschen Volkes. 



APPENDIX A: GENERALIZED FERMIONIC SUSCEPTIBILITIES 



The quantities $f . . . $^ which appear in the prefactor of the {df^n^ ■ d^n3)-teTra of the 



action Eq. ( pSj ) describing the spiral contribution to the fermion determinant are given by 
^i2{9'^9o{k-q) go,±{k-q)) , 

k 

9o,±{k-(l) 



^Yl Lio(A(A;) go,±{k)) 



J X] Lio(/^o(fc) 9o,±{k)) 
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dgo{k-q) 
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(A3) 



[9o,±{k) ^1 


1 




dgo(k~q) 


Uo,t(^) J 




^9o{k-q)\ 



9o,±{k) 



-1 



\dgo{k-qU 



dqf' 



1 9=0 



9oik) 



r%),±(fc-'?)i 



l<?=o 



|-8go,T(fc-(?) 



q=0 



9o,± 



90,T 



„2 r %)(fc~<?) i| r '9go.±(fc-9) i| 
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(A5) 



The generalized fermionic susceptibilities appearing in the ferromagnetic contribution to 
the fermion determinant (see Eq. (|108|) ) are deflned as follows: 

o 1 



9 



9 



/ - 2^o^(^) Q+{k) 



1 
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-1 



9oMk)g^^_{k) d-\k) 



-- g,\k) giik) [g' + g.^ik) g^^Lik)\ d'^k) 



(A6) 
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In the expression Eq. ( |123[ ), which is obtained for the ferromagnetic field L from the saddle- 
point solution, the following abbreviations are used: 



Ci 



Xl (Xi + 2x2 + X3) 
2(xi + X2)2(Xi + X3) ' 
X4 [(Xi - X2){Xi + Xs) + 4x1X2] 
2(xi + X2)2(Xi + X3) 
Xe , (Xt +X8)X2(X3-Xi) 



2(xi + X2) 2(xi + X2)2(xi + X3) 

Xl {xl - X1X3) X5 

(Xi + X2)^(Xi + X3) 2(xi + X3)' 



(A14) 
(A15) 
(A16) 
(A17) 
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